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Abstract 

Total cross sections for S~ and vr" on beryllium, carbon, polyethylene and copper 
as well as total cross sections for protons on beryllium and carbon have been mea- 
sured in a broad momentum range around 600GeV/c. These measurements were 
performed with a transmission technique adapted to the SELEX hyperon-beam 
experiment at Fermilab. We report on results obtained for hadron- nucleus cross 
sections and on results for (Ttot(S~N) and (Ttot(7i"~N), which were deduced from 
nuclear cross sections. 



1 Introduction 



Hadronic total cross sections provide one measure of tlie strengtli of tlie 
hadronic interaction. They have been measured for a variety of reactions over 
a broad range of center of mass energies. These studies revealed that with 
increasing center of mass (CM) energy, hadron-hadron cross sections (gener- 
ally) decrease to a minimum and then start rising again. An important current 
physics question is whether the rise of a specific hadron-hadron cross section 
is described by a power law in the CM energy. Addressing this question re- 
quires total cross-section experiments performed with a variety of hadronic 
projectiles, targets and energies covering the maximum possible range. How- 
ever, for almost 20 years, there have been few new experiments in this field. 
Thus, important hadron-hadron cross sections as o-tot(7rp) and o"tot(Kp) are 
measured only up to 380 GeV/c and the total cross section o"tot(S^p) is only 
measured up to 137GeV/c. At these maximum laboratory momenta only a 
first indication of the rise of these total cross sections is observed. 

SELEX (Fermilab E781) is a fixed-target experiment at the Fermi National 
Accelerator Laboratory using a hyperon beam of about 600 GeV/c. The SE- 
LEX spectrometer, designed for spectroscopy of charm baryons, is well-suited 
to measure total cross sections with a transmission technique. It has excellent 
scattering-angle resolution, achieved by a system of silicon microstrip detec- 
tors. 

SELEX does not have a liquid hydrogen target. Therefore, we measured 
the total hadron-nucleus cross sections o'tot(7r^Be), (Ttot(7r"C), o'tot(7r^CH2), 
o-tot(S"Be), o-tot(S"C), o-tot(S"CH2), (Ttot(pBe) and (Ttot(pC) with high preci- 
sion. We then deduced the total cross sections (Ttot(S~p) and o"tot(7i"~p) using 
both a CH2 - C subtraction technique and a method based on the Glauber 
model to derive hadron-nucleon cross sections from hadron-nucleus cross sec- 
tions. 

Further, as data on hadron-nucleus cross sections are extremely scarce for 
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charged projectiles, we also measured (Ttot(7r~Cu) and atot(S~Cu). All mea- 
surements were done during dedicated run periods in July 1997. Laboratory 
momenta range from 455 GeV/c to 635 GeV/c, the highest energy yet used for 
these studies. 



2 Experimental setup 

2.1 The hyperon beam 



The hyperon beam is generated by selecting positively or negatively charged 
secondaries around 600 GeV/c that emerge from interactions of an 800 GeV/c 
primary proton beam with a beryllium production target. Its composition has 
not been completely measured. However, we have measured the main parti- 
cle components of the event samples, which we selected to determine total 
cross sections (see section 5.2.1). This analysis shows that at the position of 
the total cross-section target the negative beam samples consist in average of 
(52.5 ± 1.6)% mesons and (47.5 ± 1.6)% baryons. Further, we measured a S~ 
fraction of (1.18 ± 0.06)% in these samples. Other baryonic fractions (p, f]^) 
were not measured, but empirical formulae (see [1]) predict they are less than 
0.1%. Likewise, the K~ fraction of the negative beam is estimated with [1] to 
be (1.6 ± 1.0)%. Thus, we expect that the vr" fraction of the event samples is 
(50.9 ± 1.9)%. 

In the event samples for positive beam we measured a meson fraction of 
(8.1 ± 1.4)% and a baryon fraction of (91.9 ± 1.4)%. Furthermore, we mea- 
sured a S^ fraction of (2.7 ± 0.7)%. Using the empirical formula given in [1], 
we expect that the tiny meson fraction consists of 70% tt"*" and 30% K"*". 

From these compositions, one sees that as long as one can distinguish mesons 
from baryons (see section 2.2), the SELEX hyperon beam offers a unique 
possibility to measure total cross sections for protons, vr", and S~ in a low 
contaminant environment. 



2.2 The section of the SELEX spectrometer used for total cross-section mea- 
surements 



The SELEX spectrometer is a 60 m long, 3 stage spectrometer. In total 
cross-section measurements, only its upstream detectors, shown in figure 1, 
are used. 

The beam spectrometer placed in front of the target, is equipped with 12 sil- 
icon microstrip detectors to track incoming particles. The first 4 microstrip 
detectors (HSDs) have a resolution (pitch/ \/l2) of 14.4 /xm and a maximum 
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Fig. 1. Sections of the SELEX spectrometer involved in the measurement of total 
cross sections. 



signal integration time of 100 ns. As this is the shortest integration time, but 
poorest spatial resolution, of all SELEX silicon microstrip detectors, the HSDs 
serve chiefly to reject stale tracks. Always, two HSDs are housed in a single 
station. The average efficiency of the HSDs is 92%. 

The remaining 8 silicon microstrip detectors of the beam spectrometer are 
grouped into 3 stations (BSSDs) mounted on a granite block inside a noise 
shielded cage (RF-cage). These detectors have a resolution of 5.8 /im and an 
average efficiency of 99.6%. 

Incoming particles are identified by a transition radiation detector (BTRD) 
with 10 separate transition radiation detector modules (TRMs). Each module 
is build of a radiator in succession with 3 proportional chambers (PCs) whose 
operating gas is a 70% Xe, 30% CO2 mixture to optimize signal response 
time and to maximize absorption of transition-radiation photons. Each cham- 
ber has a single anode readout amplifier. A radiator consists of a stack of 
200 polypropylene foils, each 17 /xm thick and spaced at 500 /xm. 
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Fig. 2. A typical BTRD signal spectrum obtained for 600 GeV/c negatively charged 
secondaries. 

Each BTRD PC gives a digital output when it detects an energy deposition 
above a fixed threshold. The sum of all PCs detecting a signal above thresh- 
old is the TRD plane count k. A typical probability spectrum of TRD plane 
counts, a BTRD signal spectrum, is shown in figure 2. It shows the baryon 
and meson responses at low and high TRD plane counts, respectively. 

The signal components are separated by fitting the function: 



Pfit(k) = E«:i(;|)pr(i-Pir'^ + E'^i(k)p^ 



■^a-Pi)'^-'^ 



baryon signal 



meson signal 



to the normalized BTRD signal spectrum. Here, pi and /tj are fit-parameters 
with the constraints 1 = /ti -|- /t2 + /ta + /t4 and pi, p2 < Ps, P4 and n is the 
maximum possible TRD plane count. The fit-parameters Pi have the meaning 
of a PC response probability, when a meson (light particle) or baryon (heavy 
particle) passes. Thus, we obtain from (1) the meson fraction (/ta -|- k/^ and 
the baryon fraction (/ti -|- k^) of the beam. 

The target is followed by the vertex spectrometer, which consists of 22 silicon 
microstrip detectors grouped into 6 stations (VSSDl, ... , VSSD5 and HSD3). 
All VSSDs have a resolution of 5.8 //m, and except for one plane, which has 
a reduced efficiency of 68%, all others have an average efficiency of 98.8%. At 
the end of the vertex spectrometer, station HSD3 is mounted to the RF cage. 

Although the total cross-section measurements presented in this article are 



based only on detectors placed in the beam and the vertex spectrometer, we 
also use other parts of the SELEX apparatus to compute corrections. Further 
detectors involved in the analysis are situated in the Ml and the M2 spec- 
trometer (see figure 1), which we describe briefiy: 

The Ml spectrometer starts at the center of the Ml magnet and ends at the 
center of the M2 magnet. For high resolution tracking of high energy particles 
in the central beam region, sets of 6 silicon microstrip detectors (LASDl and 
LASD2) are mounted to the faces of the Ml and the M2 magnet. The LASD 
detectors have a resolution of 14.4 /xm, and an average efficiency of 95.8%. For 
tracking outside the central beam region, 12 planes of wire chambers (PWCs) 
are installed. 

The M2 spectrometer starts at the center of the M2 magnet. To enhance 
the momentum resolution for high energy particles, a third station of silicon 
microstrip detectors (LASD3) is mounted to the end face of the M2 magnet. 
This station is followed by 14 PWCs that are grouped into 7 stations (M2 
PWCl, ..., M2 PWC7). 



2.3 The targets 



To optimize the precision, total cross-section measurements are done with 
special targets. Great care was taken in selecting and machining adequate tar- 
get materials in order to obtain best chemical and mechanical properties (see 
table 1). All targets are thin such that multiple scattering, quantified by ae 
of Molieres' formula is significantly lower than the 25 /irad angular resolution 
provided by the beam and vertex spectrometer. 

The carbon target is a stack of three quadratic pyrocarbon plates, each about 
5 mm thick. Pyrocarbon is composed of thin carbon layers accumulated on 
top of each other in a high-temperature methane atmosphere. Compared to 
standard graphite it offers the advantages: no open porosity, a density close to 
that of a graphite monocrystal and less than 1 ppm (parts per million) non- 
carbon constituents. The beam faces of the carbon plates were milled with a 
diamond-powder liquid and oriented such that the beam faces of the stack are 
parallel to each other. 

The polyethylene target is build from a high-purity polyethylene granulate 
with less than 1000 ppm contaminants. Molten granulate was solidified in a 
vessel, where great care was taken that no air bubbles penetrated. The ma- 
terial was then carefully machined to a target block, and beam faces were 
fiattened using a diamond pin. 

For the beryllium and the copper target, standard industry products of high 
purity are used. 



target 
material 


thickness 

L [mm] 

z-direction 


transverse 

dimensions 

X [mm] y [mm] 


density 
P* 


ere 
[/Ltrad] 


XcoU 


beryllium 


50.92 


30.7 51.2 


1.848 ± 0.002 


8.3 


16.86 


carbon 


15.46 


30.0 30.0 


2.199 ± 0.003 


6.0 


5.40 


polyethylene 


40.86 


30.0 25.0 


0.9291 ± 0.0008 


6.3 


6.66 


copper 


1.00 


30.0 30.0 


8.96 ± 0.009 


5.7 


1.05 



Table 1 

Specifics of the targets used in total cross-section measurements. L: target thick- 
ness, p*: density, ae: expected spread in scattering angle due to multiple scattering 
calculated with Molieres' formula for piab = 600GeV/c, Xcou: collision length. 

2.4 Trigger and data acquisition 



The SELEX trigger is a programmable four-stage trigger, designed to select 
events involving decays of charm hadrons in a high-intensity beam environ- 
ment. The first 3 levels: TO, Tl and T2 are hardware triggers, whereas level T3 
is an online software filter. In this section, we describe only the trigger as pro- 
grammed for total cross-section data-taking. 

At data-taking, the trigger accepted all beam events defined by the mini- 
mum-bias condition: 



TO 



51 A 52 A 53 A Kl A \/2 A V?> 



(2) 



SI, S2 and S3 are scintillation counters, and VI, V2 and V3 are veto counters 
to reject beam halo (see figure 1). In definition (2), a TO-pulse indicates a 
particle traversing the beam spectrometer in the direction of the target beam 
face. Thereby, the transverse trigger acceptance is constrained to the size of 
the rectangular hole in V2 (12.8 mm x 12.8 mm). 

In order to keep the minimum bias condition provided by the definition of 
TO, no information of detectors placed downstream of the experiment target 
infiuenced the spectrometer readout. Thus, each TO-pulse passed the Tl trig- 
ger level unbiased, and generated a T2-pulse, which started the spectrometer 
readout. The online software filter (level T3) was not used for total cross- 
section data-taking. Pulses of all trigger levels were counted by scalers for 
each spill, and saved in a trigger log file. 

The SELEX trigger controlled readout and reset of the silicon-detector sys- 
tem, the basic tool in our total cross-section measurements. Except for the 
HSDs, all other silicon detectors use an SVX-I chip technology for data readout 
[2]. SVX chips are controlled by a sequencer SRS (silicon readout sequencer) 



that interacts very closely with the trigger. First, it keeps the silicon detectors 
sensitive (for about 5/is effective integration time) and starts the chip read- 
out when receiving a T2-pulse. Second, the SRS resets the SVX-chips, when 
a silicon-clear signal arrives. Thereby, the silicon-clear signal is generated in 
the trigger logic through: 

Silicon clear = KSmuit V Cpuiscr V (Tl A T2) . (3) 



Here, Cpuiscr are pulses from a gate generator running at a frequency of 20 kHz 
and l^Smuit represents pulses generated, when the V5 veto counter (see figure 1) 
detects a high multiplicity event. The condition (Tl A T2) was irrelevant for 
total cross-section data. 



2.5 Experimental conditions and recorded data 



During the fixed-target run 1996/97, the TEVATRON was operated in 60 s 
cycles with a spill time of 20 s. Data for total cross sections were taken during 
dedicated periods, with optimized experimental conditions for this measure- 
ment. 

By adjusting the flux of the 800 GeV/c proton beam, the TO-rate was op- 
timized to run the SELEX DAQ near, but safely below its capacity limit of 
5-10^ particles per spill. The low hyperon-beam flux allowed a high silicon- 
clear rate, which resulted in a very low-noise condition for the silicon-detector 
system and a low probability for stale tracks. 

During data-taking, the Ml magnet was switched off to obtain a 2.5 m field- 
and material-free section, serving as fiducial region for precise reconstruction 
of hyperon decays. Magnet M2 was operated with a transverse momentum 
kick of p^2_ 0.84 GeV/c. 

At data-taking start, after mounting an experiment target in the RF-cage, 
an alignment RUN was taken to account for eventual detector displacements 
caused during the target installation. Then, the position of the experiment tar- 
get was alternated every 30 min between its out and in-beam position. Thus, 
almost equal amounts of data were taken with full and empty target. A RUN, 
started after each target-position change, comprised typically 10® events. A 
total of 9.8-10^ minimum-bias events were recorded with negative beam for 
the targets Be, C, Cu and CH2. With positive beam, 3.0-10^ minimum-bias 
events were written using the targets Be and C 



3 The principle of the transmission method 



In contrast to scattering experiments, where o"tot is deduced from a measured 
scattering angle distribution, in a transmission experiment o"tot is deduced from 
the number of unscattered projectiles. Strictly, unscattered means zero scat- 
tering angle, but experimental resolution and Coulomb scattering limit this to 
a determination of the number of projectiles scattered by an angle 9, which is 
smaller than a maximum angle parameter ^^ax (Fo(< ^max))- Thus, one infers 
the number of unscattered particles by extrapolating Fo(< ^max) to 6'max = 0. 

A standard transmission experiment consists of three elements: beam mon- 
itor, target, and transmission counter. The number of projectiles hitting the 
target under full-target (empty-target) condition Fq (Eq) is counted by the 
beam monitor placed in front of the target. A transmission counter, placed 
downstream of the target, counts the corresponding number of projectiles 
Ftr(< fii) (Etr(< ^i)), leaving the target within the maximum solid angles 
VLi ... r^N. Recorded counts are combined to a set of partial cross sections 
o'part(< ^i), defined as: 



0-part(< ^i) = ^-log 

pL 



Fq Et,(< l^i) 

Ftr(< a) Eo 



with p=^, (4) 



where p is the density of scattering centers in the target, A is the atomic mass 
and Na is Avogadro's number. 

Driving our choice of a transmission method is an important technical ad- 
vantage of equation (4). We do not need to know absolute efficiencies of the 
beam and the transmission monitor. Their absolute values will cancel in (4) as 
long as they remain unchanged between and during the full- and the empty- 
target RUNs (stability condition). 

Taking into account the event correlations between Fq (Eq) and Ftr(< f^i) 
(Etr(< ^i)), the statistical error of a partial cross section is given by: 
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V(<a) Fo Etr(<a) Eo 



(5) 



In a thin target approximation {pLatot <^ 1), a partial cross section apart (< ^i) 
is related to the total hadronic cross section (Jtot (see e.g. [3]) by: 



47r / > 471 



Correction for C and CN scattering 
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^i / T \ hadr i / 7 \ hadr 



elastic term inelastic term 

In equation (6), atot is infered by first correcting partial cross sections for 
Coulomb scattering (C) and the Coulomb hadronic interference (CN) and 
then extrapolating to zero solid angle. 



4 Data analysis 

4-1 Data selection 



In general, total cross-section data taken for a specific target were subject 
to varying experimental conditions: thresholds on silicon microstrip detec- 
tors, high voltages for trigger scintillators, and the inclination angle between 
primary proton beam and production target. Therefore, data belonging to a 
cross-section measurement with a specific target were divided into as many 
data sets as differing conditions had to be taken into account. This offered the 
possibility to calculate corrections and errors specifically for each experimental 
condition in a later stage of the analysis. To preserve the stability condition 
mentioned in chapter 3, a spill by spill data pre-selection was performed. Data 
of a spill or a whole run were rejected: 

(1) When the experimental conditions concerning the functionality of the 
spectrometer (detector efficiencies, trigger performance and track recon- 
struction efficiencies) suddenly changed. 

(2) When it was not possible to synchronize raw data with information in 
the trigger log file. 

(3) When the BTRD showed instabilities or when the beam phase space lay 
outside the BTRD fiducial region. 



4-2 Event selection for normalization 



The total cross-section determination is made by counting how many good 
beam tracks are removed from the beam by interactions in the target. The nor- 
malization therefore depends only on the number of good beam tracks, which 
are identified by a software decision routine. This routine reconstructs tracks in 
the beam spectrometer using the HSD and BSSD hit information. It preserves 
the minimum-bias condition for the selected data by strictly avoiding event- 
selection rules that require information from detectors placed downstream of 
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the target. An event is accepted when it is possible to reconstruct a track 
called a "norm track" , provided that the following properties are satisfied: 

(1) Not more than a total of 150 hits in all BSSDs. 

(2) At least 6 hits from BSSD planes along the track. 

(3) At least one hit from an HSD plane along the track (HSD-tagging). 

(4) A reduced track-fit x^ below 3. 

(5) An extrapolated origin of the track at the known transverse position of 
the primary production target. 

(6) Track intercept and slope parameters within the beam phase space ac- 
cepted by magnetic collimation. 

(7) A transverse track position at the longitudinal position of the experi- 
mental target, which is inside the trigger acceptance window and inside 
BTRD acceptance. 

(8) A beam momentum assigned to the track, which is ± 100 GeV/c around 
the center of gravity value of the momentum spectrum. 

Condition 3 rejects stale tracks. The selection rules 4-6 remove events in 
which hyperons decay before reaching the experiment target or react with 
detector material in the beam spectrometer. Constraint 7 assures also that 
selected tracks point to the mid-part of the experiment target face, where the 
best mechanical accuracy is obtained. 

All listed conditions were true for about 50% of the selected events. From the 
resulting set of norm tracks for full- and empty-target conditions, we establish 
classes of BTRD-tagged norm tracks. This is done by introducing cuts on the 
BTRD information as indicated in figure 2 to separate baryonic and mesonic 
norm tracks. We then determine the corresponding normalization counts Fq 
and Eq by summing the norm tracks over the appropriate signal region. 



4-3 Transmission counting 



When a norm track is found in the event, we reconstruct a single track in 
the vertex spectrometer, which is leaving the interaction target at small angle 
with respect to the norm track. 

The single-track algorithm was efficient and fast. It used hits of HSD3 to 
remove stale tracks. With loose cuts on the track parameters 98% of the norm 
tracks got assigned a track in the vertex spectrometer. Such vertex tracks were 
finally accepted as "transmitted tracks" , when: 

(1) There are at least 15 hits from VSSDs found within a track search corri- 
dor. 

(2) The reduced track-fit x^ is below 3. 
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For each transmitted track, the scattering angle 9 between norm and trans- 
mitted track is calculated. Following the idea of [4], a four-momentum transfer 
t is assigned to the event using the small angle approximation t ~ "Pbcam^^; 
where Pbeam is the momentum of the incoming particle. Transmitted tracks are 
assigned to t bins of width 5.0 ■ 10^'^ GeV^/c^. Note that we count transmitted 
tracks in t-bins, rather than in bins of solid angle Vt as discussed in chapter 3. 
Summing the events in the t-bins from zero up to a maximum t\ leads to sets 
of transmission counts Ftr(< |ti|) and Etr(< |ti|). 



4-4 Spectra of uncorrected partial cross sections 



Using the counts Fq, Eq, Ftr(< |ti|), Etr(< |ti|) and the mechanical proper- 
ties of the targets, partial cross sections o'part(< \ti\) are calculated according 




0.005 



0.015 



0.025 



0.03 



Itjl [GeV^/c^] 



Fig. 3. Spectra of uncorrected partial cross sections resulting from beryllium target 
data sets. 

to equation (4). 

Figure 3 shows some spectra for uncorrected partial cross sections. The 
strong rise of o"part(< |^i|) for \t\ < 0.002 GeV^/c^ is ascribed to multiple scat- 
tering in the target and the finite angular resolution of ~ 25/irad. Differing 
levels of partial cross-section spectra for beam particles of different kind indi- 
cate nicely the dependence of the total cross section on the projectile type. 
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4-5 Corrections for non-hadronic effects 



Partial cross sections were corrected for single Coulomb scattering (C) and 
for the Coulomb-Nuclear interference effect (CN) evaluating the expression: 

<:;« i*>i) = .^,« H) -/ (I) *' -/ (!) ^*' . (7) 



C correction CN correction 

Applying the Coulomb correction, a change in the extrapolated cross section 
of not more than 0.5% is observed for the light targets Be, C and CH2. For the 
Cu target a change of up to 11% is noticed. The CN correction is roughly one 
order of magnitude smaller than the Coulomb correction and has negligible 
effect on the extrapolated cross section. 



J^.6 The extrapolation method 



As |t| approaches zero, the growth behavior of partial cross sections is gov- 
erned ideally by the elastic term in equation (6). At small \t\, hadronic coherent 
elastic scattering off nuclei dominates. Thus, we obtain for the elastic term in 
equation (6), the expression: 



( ( da\ 



hadr 



CT, 



dt' = —^^(1 + pS 



rft'y^i iq-rB, 



nuc 



1 - e 



-Dnuc^^ 



where -Bnuc is the exponential slope observed in hadronic coherent elastic scat- 
tering off nuclei. Therefore, we choose the functional form 



/(ai,a2,t) = ai 1 - e"^' (9) 



to describe the variation of partial cross sections with respect to |ti|. 

The parameters cti and 0:2 are determined in fitting function (9) to differ- 
ences in corrected partial cross sections of adjacent t-bins. Thereby, only those 
corrected partial cross sections where ti is in the range of tmin = -0.007 GeV^/c^ 
to tjnax = -0.03 GeV^/c^ enter the fit procedure. The limits tmax and t-^in ac- 
count for experimental sensitivity to hadronic coherent elastic scattering off 
nuclei. Their derivation is described in section 4.6.1. 

Starting from the partial cross section apart (< l^minl), the extrapolation to 
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the total cross section o"tot is determined by accounting for the expected growth 
in partial cross sections from tmin to t = using the expression: 



C"tot — Crpart(< KminI) + Ctl 



1 - e' 



a2t 



min 



(10) 



4-6.1 The limits tmin and t^^x o-nd the sensitivity of the SELEX experiment 
to coherent hadronic elastic scattering off nuclei 

In measurements of hadron-nucleus cross sections, it is essential that the ex- 
periment is sensitive to hadronic coherent elastic scattering off nuclei. Further, 
one must be able to distinguish coherent from incoherent scattering processes 
off nucleons. In scattering off nuclei, the nucleus can break up when the energy 
transfer exceeds the binding energy of its nucleons. This leads to a contribu- 
tion of incoherent scattering off nucleons for |t| > 0.015 GeV^/c^. In that 
case, the hadronic differential elastic cross section, entering the elastic term 
of equation (6), contains two parts: 

\ dt , IGtt IGtt 

\ / ci V ^ -■ V ::,^ ^ 

coherent scattering incoherent scattering 

There is a term for coherent elastic scattering off the nucleus, in which o"tot(hA) 
is the total nuclear cross section, and a term for incoherent scattering off nu- 
cleons, in which o"tot(hN) is the corresponding hadron-nucleon cross section. 
i?N is the slope parameter for scattering off nucleons, and N(A) is a factor de- 
scribing the effective number of nucleons taking part in the incoherent process 
for target nuclei of mass A (see [5]). 

The contribution of the incoherent term decreases the growth behavior of 
the elastic term in (6) because B^ is typically one or more orders of magnitude 
smaller than -Bnuc- An extrapolation based on partial cross sections, selected 
in a |t|-range far above 0.015 GeV^/c^ would lead to a systematically lowered 
cross-section result as a fraction of the elastic processes is not discriminated. 
Consequently, we looked for a t-interval [tmin; ^max] to select partial cross sec- 
tions where their growth is dominated by -Bnuc- 

The sensitivity of the SELEX spectrometer to hadronic coherent elastic 
scattering off nuclei was verified by looking at background subtracted but not 
acceptance corrected differential scattering spectra, which are defined by: 



S(t) 



pLT 



Fn En 



(12) 



Here, F is the width of the t-bins. F(t) (E(t)) is the number of scattering 
events found in the full-target (empty-target) data sets that fall into the in- 
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terval [|t| -T/2 ; \t\ +T/2]. 

Figure 4 shows a typical example of an S(t) spectrum obtained for S~ scat- 
tering off carbon nuclei. The spectrum shows three regions governed by ap- 
parently different exponential slopes, which can be explained by contributions 



^ 10 




0.35 
t [GeV^/c^] 



Fig. 4. Differential scattering spectrum obtained for S carbon reactions, showing 
the Coulomb, the coherent and the incoherent region. 

of Coulomb scattering, coherent elastic scattering and incoherent elastic scat- 
tering comparable to measurements described in [5]. 

Determinations of the slope parameters -Bnuc and Sn in S(t) spectra showed 
the expected order of magnitude for all targets, and -Bnuc agreed quite well with 
data presented in [6] . Furthermore, the magnitude of i?nuc is also reflected by 
the size of parameter 0:2 in equation (10), when applying the extrapolation. 

From such studies, we choose tmax = -0.03 GeV^/c^ fixed, as this value is 
well inside the region dominated by coherent hadronic elastic scattering off 
nuclei for all targets. The contribution of the integrated incoherent term at 
this tmax is much lower than the integrated coherent term. 

Further, to avoid large multiple-scattering corrections, we choose a tmin of 
—0.007 GeV^/c^, such that the angular resolution has negligible effect on the 
extrapolated total cross section. 
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5 Corrections 



5.1 Trigger-rate corrections 



The trigger rate influences tlie reconstruction efliciency for tracks and thus 
alters the transmission ratios Tfuu and Tempty per spilL Figure 5 shows an 
instructive example of this effect. 

Due to the rate effect, our extrapolated total cross-section experiences a 
shift Ato when the average TO-counts, calculated for all empty and all full- 
target spills separately, differ. 
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Fig. 5. Dependency of full-target transmission ratios on the TO-count. 

To determine the shift Ato we calculate full and empty-target transmission 
ratios per spill for |t| < 0.01 GeV^/c^ and describe their rate dependency by 
fitting to the expression 



Tfit(TO) = ;3i,k + ,32,kT0'^ . 



(13) 



We have studied the effect of different powers (k = 2, 3, 4) to estimate sys- 
tematic errors. 

We choose the average TO-rate TO, comprising all full and all empty-target 
spills as reference rate for the rate correction. Thus, transmission ratios per 
spill are corrected by evaluating: 



corrected 



Tj(|t| < O.OlGeVVc^) + /32,k(T0 -TOf 



uncorrected 



(14) 



correction 
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which results in a set of corrected transmission ratios T™. Fit-function de- 
pendent offsets Axck are deduced by: 

ATo,k = a^'i{< O.OlGeVVc^) - apart(< O.OlGeVVc^) (15) 

and averaged to a mean offset Ato- Total cross sections are then corrected by: 

c^to? = ^tot + Ato • (16) 

Averaged sizes of the rate correction are presented in table 2. We want to 
mention that the copper data were taken at higher rate, where the slope of 
function (13) is steeper. This, together with the small thickness of the copper 
target, causes large corrections. 

5.2 Corrections for beam contaminants 



A transition radiation detector does not make an exact particle identifica- 
tion because of statistical fiuctuations in X-ray generation and background 
from various processes. Therefore, when selecting the baryon or the meson 
component of the hyperon beam by applying cuts on the BTRD plane count, 
we need to account for: 

(1) The amount of meson (baryon) contaminants in the baryon (meson) sam- 
ple and the effect on the total cross section. 

(2) The amount of baryon (meson) contaminants in a specific sample for a 
measurement with protons or S~ (tt") and the effect on the total cross 
section. 

Once the contaminant fraction e is determined, the experimental cross section 
(Tt^f can be corrected by the term Acont using: 



pL 



A2)(^-pL(aill -411) 



(17) 



Correction Ac 



This formula was derived in [7] for a two component beam having a contami- 
nation fraction e*^^-*. 



5.2.1 Beam contaminant determination 

In a first step, total cross sections resulting from data sets are corrected 
for the fraction of mesons (baryons) in a baryon sample (meson sample) us- 



ing (17). Therefore, we fit function (1) to normalized BTRD signal spectra, 
which are recorded for norm tracks. For negative beam these fits yield an 
average baryon fraction (/ti + K2) of (47.5 ± 1.6)% and an average meson 
fraction (^3 + K4) of (52.5 ± 1.6)%. For positive beam, we measure a baryon 
fraction of (91.9 ± 1.4)% and a meson fraction of (8.1 ± 1.4)%. To deduce 
the meson (baryon) contaminant fraction e, we sum the meson (baryon) com- 
ponent of (1) over the TRD plane count region shown in figure 2. Further, 
the difference al^l — alol is calculated in taking rate corrected extrapolated 
cross-section results obtained for the meson and the baryon beam component. 
In a second step, we account for the main contaminant disturbing a specific 
measurement for protons, S~ and vr^. According to the expected hyperon- 
beam composition we correct: 



(1) For the effect of H~ particles in the baryon sample, when measuring S^A 
cross sections. 

(2) For the effect of S+ particles in the baryon sample, when measuring pA 
cross sections. 

(3) For the effect of K~ particles in the meson sample, when measuring vr^A 
cross sections. 

For case (1), we measure the overall fraction of S~ particles in each negative- 
beam data sample and for case (2) we measure the overall fraction of E+ parti- 
cles in each positive-beam data sample. Therefore, we count the decays E~ — >■ 
n + 7r~, S~ — »• A° -|- 7r~ and S"*" — > n -|- 7r+, reconstructed for a known amount 
of norm tracks within the field-free region of the Ml magnet. Figure 6 shows 
some hyperon-mass spectra obtained by the decay reconstruction. 
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Fig. 6. Hyperon-mass spectra obtained from reconstructed E~, H~ and S^ decays. 
The spectra are fit to a Gaussian plus a linear background function, mx is the mean 
mass found for hyperon X and (Tm is the corresponding mass resolution resulting 
from the fit. 
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Particle decay counts are corrected for geometrical acceptance, branching ratio 
and decay losses after the target, to yield the overall hyperon contaminant 
fractions. Here, we find an overall S~ fraction of (1.18 ± 0.06)%, and an 
overall S+ fraction of (2.7 ± 0.7)%. These fractions are then divided by the 
baryon fraction (/ti + k,2), known from the first step procedure to yield the 
hyperon contaminant fraction e of the baryon component. 

Case (3) requires knowledge of the number of K~ particles in the meson 
sample. As the SELEX spectrometer does not provide a tool to discriminate 
600 GeV/c 7r~ against K~ particles, we estimate the overall fraction of K~ 
particles in the sample using particle-flux parameterizations of [1]. This results 
in an overall K~ fraction of (1.6 ± 1.0)%, which divided by the meson fraction 
(k3 + K4) yields the K~ contaminant fraction of the meson component. 

Calculating the contaminant correction using equation (17) requires know- 
ledge of the total cross sections (Ttot(S~A), crtot(S^A) and (Ttot(K A). As data 
on these cross sections are either scarce or do not exist, we estimate them 
using approximations like: 

,„_ . ,^_ fftot(pA) 

o-tot(^ A) ^ o-tot(^ p) 7 — r , (18) 

^tot(pp) 



and neglect weak but existing energy dependencies. Necessary data for hadron- 
nucleon cross sections are taken from [4,8] and data for pA-cross sections are 
taken from [9]. 

Averaged sizes of the contaminant correction including both correction steps 
are shown in table 2. 



6 Results for hadron-nucleus cross sections 



Total cross sections as well as their statistical and systematic errors were 
determined for each dataset separately. In order to calculate average total cross 
sections and average systematic errors, we use weighted means. We present 
the error contributions, the data averaging method and the final results. 



6.1 Measurement errors 



6.1.1 The statistical error 

The dominant error contribution is the statistical error, which is governed 
by the statistical uncertainty of the partial cross section o"part(< l^minl), used in 
the extrapolation. Further statistical error contributions, originating in other 
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terms of the error propagated formula (10), are negligible. The statistical errors 
for each measurement are presented in table 3. 



6.1.2 Systematic errors 

In this section we briefly describe the systematic errors found during the 
data analysis. Table 2 gives an overview of the average sizes of these errors as 
well as the rate correction and the contaminant correction. 

Systematic error of the extrapolation 5*^^*'' 

A significant error contribution is the systematic error in the extrapolation 
of partial cross sections. This error takes into account the RMS-spread (root 
mean square) in the extrapolated total cross sections with respect to the ex- 
trapolated total cross section at tmin = — 0.007 GeV^/c^, when tmm is varied 
from -0.004 GeVVc^ to -0.01 GeVVc^. 

Cut on the BTRD signal spectrum ^^^rd 

Although contaminant and rate corrections are applied for each specific cut 
on the BTRD signal spectrum, we still observe a variation of the cross section 
when varying the cut on the TRD plane count by ± 1 unit around its nomi- 
nal value. Therefore, we calculate a systematic error, which is the maximum 
spread in the cross sections found in the cut variation. 

Spill to spill fluctuations 5^^'^ 

Here, we compare the statistical error in crpart(< 0.01 GeV^/c^), which we 
calculate by (5) with the error in o"part(< 0.01 GeV^/c^) calculated from the 
experimentally observed RMS-spread of rate corrected transmission ratios per 
spill. The difference in these errors accounts for remaining non statistical spill 
to spill fluctuations. 

Systematic error of the rate correction 5^^^*^ 

This error takes into account the error arising from different functional at- 
tempts to describe the rate effect presented in section 5.1. Its value is given 
by the maximum spread of the ATo,k with respect to their average value Ato- 

Systematic error of the contaminant correction 5'^°°* 

This systematic error accounts for the uncertainty in the fit parameters of 
the four-fold binomial distribution (1) and for the uncertainty in the contam- 
inant fractions for S+, S~ and K~. 

Uncertainty of the target density 5*^' 
The target densities were measured several times, using a pycnometer and 
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a buoyancy method. Laboratory studies showed systematic discrepancies in 
the density measurement, which are included in the density errors shown in 
table 1. These errors are propagated to an error contribution to the total cross 
sections, which are on a 0.1% level. 



cross 

section 


Plab 

[GeV/c] 


^cxtr 

[mb] 


^BTRD 
[mb] 


stemati 

[mb] 


c errors 

^ratc 

[mb] 


^cont 

[mb] 


5tgt 

[mb] 


Corrections 

Ato Acont 

[mb] [mb] 


CTtot (pBe) 


536 


0.91 


0.70 


0.25 


0.35 


0.06 


0.30 


-1.24 


0.62 


o-tot(S"Be) 


638 


1.20 


0.49 


0.04 


0.10 


0.07 


0.27 


-0.93 


0.65 


o-tot{7r"Bc) 


638 


0.50 


0.17 


0.22 


0.05 


0.61 


0.21 


-0.79 


1.00 


crtot (pC) 


457 


0.90 


2.11 


0.54 


0.38 


0.09 


0.47 


11.22 


0.91 


CTtot (pC) 


490 


1.81 


1.53 


0.68 


1.15 


0.10 


0.47 


-3.87 


0.86 


atot(S-C) 


598 


1.57 


1.92 


1.21 


1.18 


0.13 


0.43 


-6.42 


1.12 


CTtot (tt" C) 


591 


1.30 


1.40 


1.50 


0.95 


0.63 


0.33 


-3.11 


1.03 


(Ttot(S-CH2) 


589 


2.10 


2.55 


0.69 


0.16 


0.16 


0.30 


3.67 


1.44 


o-tot(7r"CH2) 


585 


1.26 


0.96 


0.54 


0.12 


0.75 


0.23 


2.90 


1.21 


^tot(S-Cu) 


609 


163 


41 


76 


41 


0.33 


1.23 


-754 


3.1 


o-tot(7r"Cu) 


608 


85 


52 


78 


36 


2.99 


1.03 


-649 


4.7 



Table 2 

Average sizes of systematic errors and corrections. For explanation of symbols see 

text of sections 5.1, 5.2 and 6.1. 



6.2 Data-averaging and results on hadron-nucleus cross sections 



6.2.1 The average total cross section 

Total cross-section results (Jtot,i, obtained from i = 1, ..., N data sets, are 
combined to an average total cross section atot with a statistical error ^'^^'^^o'tot 
and an average systematic error 5^^^^ a tot. The results are shown in table 3. 

We average the total cross-sections (Xtot.i that correspond to a specific mea- 
surement using the weighted mean: 



C^tot 



E ^i <7tot,i 

i=l 

N 

E^i 

i=l 



0J\ 



{61 



Stat \ 2 



M 

j=i 



syst\2 



(19) 



The weight uo^ includes the statistical error ((5^*^^*) of data set i and all j = 1, 

..., M systematic errors S^i" described in section 6.1.1 and 6.1.2. 

As the statistical error is supposed to decrease, when adding more data to the 
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evaluation, we calculate the statistical error in the averaged total cross section 
by: 



^stat^ 



tot 



1 



MVE(^stat)2 • 



(20) 



In assigning a systematic error to an average total cross section atot we assume 
that the systematic errors of the single measurements can be just averaged. 
Therefore, we quote as average systematic error: 



^syst 



Ctot 
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[m ^ 1 


E^i 


E(0' 
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Lj=i J 



N 
i=l 



(21) 



which accounts for the weights uj\ used in (19). 

Further, we quote a total error (5*°*atot of the average total cross section, 
which is calculated from: 



<5*°*^tot = \/(5^*^*^tot)' + (5^^^*^tot) 



f22) 



cross 
section 


Plab 

[GeV/c] 


O"tot 

[mb] 


<5^'^'atot 
[mb] 


<5^^^*^tot 

[mb] 


[mb] 


crtot(pBe) 


536 


268.6 


±0.7 


± 1.3 


± 1.5 


atot (S- Be) 


638 


249.1 


±0.9 


± 1.3 


± 1.6 


crtot(vr"Be) 


638 


188.7 


±0.8 


±0.9 


± 1.2 


o-tot(pC) 


457 


333.6 


±3.1 


± 2.4 


± 3.9 


o-tot(pC) 


490 


335.4 


±3.6 


± 2.9 


±4.6 


^tot(S-C) 


598 


308.9 


± 2.1 


± 3.8 


±4.3 


o-tot(7r^C) 


591 


234.1 


± 1.5 


± 3.1 


± 3.5 


C7tot(S-CH2) 


589 


376.4 


± 2.0 


± 4.1 


± 4.5 


o"tot(7r~CH2) 


585 


286.1 


± 1.3 


± 2.0 


± 2.4 


atot(S-Cu) 


609 


1232 


± 133 


± 192 


± 233 


o-tot(vr"Cu) 


608 


1032 


± 77 


± 162 


± 179 



Table 3 

Results for nuclear total cross sections. For explanation of symbols see text of sec- 
tion 6.2. 
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6.3 Comparison to existing data on hadron-nucleus total cross sections 



A literature survey showed that experimental data on hadron-nucleus to- 
tal cross sections for charged projectiles at high energies is extremely scarce. 
Information for proton-nucleus and 7r~-nucleus total cross sections is only pro- 
vided by the references [5,6,9] and displayed together with our results in the 
figures 11, 7, 8, 9 and 10. No data were found for S~-nucleus total cross 
sections. 



6.3.1 Comparison of nucleon-nucleus total cross sections 

In the figures 7 and 8 we display a compilation of proton-nucleus and 
neutron-nucleus cross sections extracted from [7,10-12] together with our re- 
sults. As can be seen, the proton-nucleus cross sections of [5] at piab = 20 GeV/c 
and the neutron-nucleus cross sections are close by. For this reason we assume 
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Fig. 7. Summary of experiment data on (Ttot(nBe) from [7,10-12] and on crtot(pBe) 
from [5] and SELEX. Overlaid are results from the model calculation (see section 7). 

that differences between neutron-nucleus and proton-nucleus cross sections 
are negligibly small above 20 GeV/c. This allows a comparison of our proton- 
nucleus cross sections with corresponding neutron-nucleus cross-section data 
available at much higher energy. 
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Comparing our results with neutron-nucleus cross sections at 131 - 273 GeV/c 
(data of [7]) shows that our measurements join the sequence of these data 
points. Averaging the neutron-beryllium total cross sections in this momentum 
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Fig. 8. Summary of experiment data on crtot(nC) from [7,10-12] and on o"tot(pC) 
from [5] and SELEX. Overlaid are results from the model calculation (see section 7). 

range results in 271.0 ± 0.6mbarn, which is close to our proton-beryllium 
cross section at 536 GeV/c of 268.6 ± l.Smbarn. A similar calculation for 
the neutron-carbon cross section gives a mean value of 331.0 ± O.Smbarn, 
which is close to our measurements of the proton-carbon cross section around 
457GeV/c of 333.6 ± 3.9mbarn. 



6.4 Comparison of n -nucleus total cross sections 



High-energy data for atot(7r^Be), atot(7r~C) and atot(7r^Cu) that were deter- 
mined using a transmission technique are presented in the thesis of A. Schiz [9]. 
Unfortunately, the statistical errors quoted for the 7r"A total cross-sections are 
too large and miss further corrections. Therefore, we present in figure 11, 9 
and 10 "normalization factors" for 7r~-nucleus scattering taken from [6]. These 
factors are based on the thesis [9], have the meaning of a total 7r~-nucleus 
cross section, but result from fits to 7r~-nucleus elastic scattering data, and 
have smaller error. 
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Fig. 9. Summary of experiment data on (Ttot(7r Be) from [6] and SELEX. Overlaid 
are results from the model calculation (see section 7). 
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Fig. 10. Summary of experiment data on (Ttot(7r C) from [6] and SELEX. Overlaid 
are results from the model calculation (see section 7). 
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Fig. 11. Results for CTtot(vr~Cu) from [6] and SELEX. 

The figures 9, 10 and 11 show that the SELEX resuhs for 7r~Be, 7r~C and 
7r~Cu cross sections join the displayed normahzation factors. However, further 
and more precise data would allow a more detailed comparison. 



7 Model description of hadron-nucleus cross sections 



In this section, we introduce a model calculation for hadron-nucleus cross 
sections and show how well it describes the data. 



7.0.1 The Glauber model and the inelastic screening correction 

As shown in [7], the Glauber model [13,14] including an inelastic screening 
correction [15], is very precise in describing neutron- nucleus cross sections at 
high energy. The Glauber model accounts for the elastic screening effect in 
nuclei via multiple elastic scattering between the incident hadron h and the 
nucleons N. As mentioned in [7], nuclear total cross sections calculated by 
the Glauber model exceed experimental data. This is compensated by taking 
into account the inelastic screening correction described in [15]. It accounts 
for inelastic reactions h + N —^ N + X, which produce an inelastic screening 
effect. Consequently, a model cross section 0"™°'^ comprises two parts: 



<f(A,atot(hN)) 



a. 



Gla 



(A,atot(hN)) 



Act 



Kar 



(23) 



These are a Glauber model cross section ^{^^.^(A, o"tot(hN)) and an inelastic 
screening correction Acr^^'". 
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The Glauber model cross section 
According to [14], o-^J.'^(A, o"tot(hN)) can be calculated by: 



1 - ^^^^a,o,(hN)T(b)'^ 



bdb } (24) 



{cxj 
/l - 

oo oo 

T(b) = — / Jo(qb)e-^™^g('q)qjq ^ gj'q) = — / r sin (qr)p(r)dr . 
27r J q J 

o o 

Here p'^^ i^ ^^e real to the imaginary part of the elastic scattering amplitude 
in the forward direction observed in hadron-nucleon elastic scattering and b 
is the impact parameter. i?hN is the hadronic slope in hadron-nucleon elastic 
scattering and Jo is a Bessel function of O'th order. The nuclear density p(r) 
is normalized as: 

oo 

47r f p^iydr = 1 . (25) 



The inelastic screening correction 

The inelastic screening correction Aa^""^, originally formulated in [15] for 
proton-nucleus reactions, is generalized by: 



oo {v^-mp)2 

d^cr 



Aa^- = 47r / / ^^ g-i.,.,(hN)At(b) p(q^^ g) ^ ^^2^2^ (26) 
J J VdtdM^/^ „ 

+00 +00 

t(b)= I p{h,z)dz , F(qL,b) = A / p(b, z)e^'iL-dz 



—00 —00 

qL = (M^ - mj)^ . 
^ s 



Here nip is the proton mass and m.,^ is the pion mass. The double differential 
cross section d^a/dtdM^ describes the inelastic reaction h + N^N + Xoi 
the incident hadron h with a nucleon A^, where the resulting final state X has 
an invariant mass squared of M^. 



7.0.2 Input parameters for the total cross-section model 

Model input parameters are crtot(hN), p'^^^, -BhN, />(!■) and (d^(T/dtdM^)|t=o- 
All of them are extracted from experimental data with A^ = p. 
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Input parameter o"tot(hN) 

To be able to access o"tot(hN) for special ranges of center of mass energies 
y^, data on pp and Tr^p total cross sections, taken from [8], were fit to a 
smooth function: 



crtot(hp,s) = ^ + a2log^(s) 



(27) 



The fit-parameters Oi, their errors and the validity range of each parameteri- 
zation, are shown in table 4. The result of each parameterization is in mbarn, 
when using s in GeV^. 



reaction 


ao 


ai 


&2 


momentuin range 


PP 

vr^p 


49.51 ± 0.26 
55.2 ± 7.2 


0.097 ± 0.002 
0.255 ± 0.032 


0.314 ± 0.004 
0.346 ± 0.020 


10 ... 3000 GeV/c 
80 ... 380 GeV/c 



Table 4 

Fit-Parameters and validity range of the total cross-section parameterizations. 



Input parameter p(^p 

We parameterize Ppp(piab) and p^_p(piab), using data on Ppp from [16-28] 
and data on p^_ from [16,21,29,30], assuming that p' reaches a constant value 




10 10 

Piab [GeV/c] 



P,ab [GeV/c] 



Fig. 12. Our parameterizations for Ppp(piab) and p^-p(piab) together with experi- 
mental data from [16-28]. 
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when piab goes to infinity. Our fits are 



Ppp(P 



6i 



labi 



„0.742 
Plab 



6.6 

37599 
Plab 



0.124 



for 0.8 GeV/c < piab < 2100 GeV/c 



p'tt-pIp: 



lab J 



0.92 

'0.54 
Plab 



+ 0.54 



(2J 



(29) 



for 8.0 GeV/c < piab < 345 GeV/c , 



wliere piab is in GeV/c. Figure 12 displays tliese fit-functions togetlier witli all 
data points included in the fit. 



Input parameter B] 



hp 



For the hadronic slope parameters -Bpp and -Bjr-p we take the parameteriza- 
tions presented in [31]: 



Bpp(p 



pp I Plab J 



Bpp,i 

Bpp,2 
Bpp,3 



11.13 
9.26 
9.67 



1.21 



/Plab 



0.301og{piab} 
^ + 0.28 log{piab} 

^ + 0.101og{piab} 



q2 = 0.02 

q2 = 0.20 (30) 

q2 = 0.40 



B,rp(P 



lab J 



Bjrp,! 
B7rp,2 
Bjrp.S 



9.11 + 

6.95 + 
6.13 + 



0.65 

VPlab 

0.65 

VPlab 

0.65 

•\/Plab 



+ 0.291og{piab} 
+ 0.271og{piab} 
+ 0.251og{piab} 



q2 = 0.02 

q2 = 0.20 (31) 

q2 = 0.40 . 



Here, q^ is in units of GeV^/c^. These parameterizations are linearly interpo- 
lated to account for the dependency of BhN on both piab and q^. 



Input parameter (d^cr/dtdM^ 



|t=o 



To calculate the inelastic screening correction Acr , we use the parameter- 
ization of (d^o"/dtdM^)|t=o for the process p + p ^ p + X, given in [7]: 



UtdM2A=o 



'26.470(M2-1.17) - 35.969(M2 - 1.17)^ + 
18.470(M2 - 1.17)=^ - 4.143(M2 - 1.17)^ + 
0.341(M2 - 1.17)5 for 1.17 < M^ < 5 GeVVc^ (^2) 

for M2 > 5 GeVVc^ 



4.4 
M2 



In addition, we also use more recent parameterizations for (d^(T/dtdM^)|t=o 
to describe the processes p + p ^ p + X and tt + p ^ p + X, which are pre- 
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sented in [32] and are based on calculations of triple- Regge diagrams in [33]. 
For M^ < Mq, these parameterizations consist of a background term and a 
sum of non-energy-dependent resonance terms. In case M^ > M^ the param- 
eterizations consist of a sum over contributions from triple- Regge diagrams: 



' dV 

dtdM2 



t=o 



4^ (M2-Mf)2+ri "^ (M2-M;^;^)2+df 



S^~"k(0) 



M2 <M2 
M2 >M2 



(33) 



Instead of displaying the large amount of parameters for equation (33), which 
are taken from calculations in [33], we display the parameterizations (32) 
and (33) in figure 13. 



^7.5 



-r 5 









2.5 



p + p -> p + X 
p =600 GeV/c 



'lab 




Function (32) 

Function (33) 
for pp reaction 



15 



20 



M^ [GeV^/c"*] 



."6 

> 

« 
14 






7t' + p ^ P + X 

Piab=600 GeV/c 

Function (33) 
for Tip reaction 



Mf 



10 



15 



20 



M^ [GeV^/c"*] 



Fig. 13. The parameterizations (32) and (33) evaluated for piab = 600GeV/c. 

Compared to (33), parameterization (32) has no s-dependence. Further, pa- 
rameterization (33) is not continuous and the resonance sizes are quite different 
for p + p — i> p + X. 

Input parameter p(r) 



In the calculations, we use density distributions p(r) that are based on the 
harmonic-oscillator model: 



p(r) 



Po 



1 + al^-V 



^rad ' 



(34) 



This offers the possibility to calculate some integrals in an analytic way and 
gives a better description of the (charge-) density distribution for light nuclei 
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than a standard two-parameter Fermi parameterization. 



nucleus 



data from [34] 
e-A scattering 
arad [fm] 



a 



fit result using (32) 

ina-°''(NA,arad) 
arad [fm] 



fit result using (33) 

in<?''(NA,arad) 

arad [fm] 



beryllium 
carbon 



0.611 
1.067 



1.791 
1.687 



1.89981 
1.79247 



2.02914 

1.89277 



Table 5 

Parameters of the density distribution /5(r) from electron-nucleus elastic scatter- 
ing [34] and the radius parameters resulting from a fit of (T™°'^(arad5Piab) to nA 
cross-section data in the momentum range 10 - 273GeV/c. 



As reported in [7], we also find that the model does not provide a good 
description of neutron-nucleus total cross sections if one uses both a and 
arad from electron-scattering data [34]. Therefore, we used a values from [34] 
and adjusted the radius parameter arad, such that the model cross section 
cr™°'^(NA, arad) gives a best description of nA-cross section data in the mo- 
mentum range 10 - 273 GeV/c. Adjusting of arad was done for each nucleus 
and for each of the parameterizations (32) and (33) separately. Table 5 gives 
a summary of the density parameters. 



1.0.3 Results of the model calculations 

Results for nucleon-nucleus model cross sections 

To show the quality of our model calculation after adjusting the nuclear den- 
sity parameter arad, we evaluated the total cross sections o"™°'^(Be, o"tot(pp)) 
and o-™°'^(C, o"tot(pp)) using function (32). This was done for data on o"tot(pp) 
taken from [8] and for values on (7tot(pp) resulting from our fit (27). The cal- 
culations were done at many different values of piab to show the behavior over 
the entire high momentum region. Scatter in the model calculations (observed 
when experimental data on (Ttot(pp) are used) demonstrate the sensitivity of 
the model to small changes in o"tot(pp)- 

Summaries of calculation and data are shown in figure 7 and 8. They show 
that the calculations refiect quite well the cross-section data for piab > 5 GeV/c. 
The nBe data of [7] in the range 131 - 273 GeV/c suggest a rise of the nBe 
cross section with energy that is also indicated by the model calculation. Our 
data point does not show any rise for pBe. In the case of nC cross sections our 
measurements join both data at lower energy and calculation very nicely. 

Results for 7r~-nucleus model cross sections 

We evaluated the cross sections o"™°'^(Be, o"tot(7r^p)) and (tJ^°'^(C, (Ttot(7r^p)) 
using function (33) and the corresponding nuclear density parameter arad, 
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which was determined by a fit of the model cross section to neutron-nucleus 
data. All further input parameters are specific for 7r~p-reactions. The calcu- 
lations were done for data on 0"tot(7i"~p) taken from [8] and for values from 
function (27). 

Results are shown in figure 9 and 10 together with data for 7r^-nucleus total 
cross sections from [6] and the SELEX experiment. The figures show that the 
calculations match our measurements quite well and agree within errors with 
lower-energy data from [6]. 



8 Results for hadron-nucleon cross sections 



The hadron-nucleon cross sections o"tot(S~N) and o"tot(7r~N) were first deter- 
mined by a CH2 - C method. As this method provides hadron-nucleon cross 
sections only with a precision on the order of 10%, we improved the precision 
using a method which takes advantage of the more precise hadron-nucleus 
cross-section ratios. 



?. i Hadron-nucleon cross sections using a CH2 - C difference method 



The hadron-nucleon cross sections o-tot(S N) and o"tot(7r N) can be deduced 
from corresponding cross sections measured on carbon and polyethylene by: 

atot(hN) = i[(Ttot(hCH2) - atot(hC)] , (35) 



where h denotes the incident hadron. Results obtained by this method are 
presented in table 9. The quoted errors are calculated from the total errors in 
the hadron-nucleus cross sections given in table 3. 



3.2 Hadron-nucleon cross sections deduced from hadron-nucleus cross sec- 
tions 



In a second approach, we deduce hadron-nucleon cross sections from ratios 
of measured hadron-nucleus cross sections. To motivate the method, we first 
derive empirical relations between hadron-nucleon and hadron-nucleus cross 
section ratios, which we then refine using the model calculation described in 
section 7. 

To derive empirical relations between hadron-nucleon and hadron-nucleus 
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cross section ratios we use data on hadron-nucleon cross sections around 
137GeV/c from [4,8], and obtain the hadron-nucleon cross-section ratios: 



o-tot(pp) 



0.635 ±0.006, 



Q-totlS p) 

(^tot(PP) 



0.901 ±0.012 



(36) 



Next, we build nuclear cross-section ratios using our measurements for the 
S^A, TT^A and pA cross sections from table 3. 

Our pA cross sections were measured at lower laboratory momentum than 
the corresponding S^A or vr^A cross sections. To correct for this, we scale the 
pA cross sections by a factor kgcaie before building the cross-section ratio. The 
scale factor takes into account the growth of the pA cross section from the 
laboratory momentum where it was measured to the larger laboratory mo- 
mentum of the corresponding S^A or vr^A cross section. Scaling factors are 
calculated using the model described in section 7. They are displayed together 
with the nuclear cross-section ratios in table 6. 



scaled cross-section ratio 


Piab [GeV/c] 


J^scalc 


atot(^~Be)/cJtot(pBe) = 0.698 ± 0.006 


640 


1.0058 


cTtot(vr-C)Mot(pC) = 0.695 ± 0.014 


590 


1.0036 


c7tot(S-Be)/c7tot(pBe) = 0.922 ± 0.008 


640 


1.0058 


c7tot(S-C)Mot(pC) = 0.917 ± 0.018 


590 


1.0040 



Table 6 

Nuclear cross-section ratios. The pA-cross section is scaled by kgcaic to account for 

the discrepancy in laboratory momenta of the cross sections used in the ratio. 

The nuclear ratios show that the 7r~A cross sections are about 0.7 times and 
the S^A cross sections are about 0.92 times smaller than the pA cross section. 
To get a first relation between hadron-nucleon and hadron-nucleus cross 
sections, we ignore the weak energy dependence of the cross-section ratios. 
Calculating the ratios of hadron-nucleon to hadron-nucleus cross-section ratios 
using the above data gives the results presented in table 7. 



The double ratios show a small but significant deviation from one especially for 
ratios involving 7r~ cross sections. From this empirical observation it follows 
that a hadron-nucleon cross section atot{hN) can be approximately derived 
from the pp cross section and a hadron-nucleus cross-section ratio using the 
relation: 



(Ttot{hN) ^ K X CTtot(PP) X 



Vtot(hA)^ 
.o-tot(pA)^ 



(37) 
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double ratio 


result 


double ratio 


result 


CTtot (S^ p)/o-tot (pp) 

CTtot(S-Bc)/atot(pBc) 

crtot(S"p)/o-tot(pp) 
CTtot(S-C)/(Ttot(pC) 


0.977 ± 0.016 
0.983 ± 0.023 


CTtot (7r^p)/o-tot (pp) 
o-tot (tt- Bc)/o-tot (pBc) 

o-tot(7r"p)/o-tot(pp) 

CTtot (tT^ C)/o-tot (pC) 


0.910 ± 0.012 
0.915 ± 0.020 


average (k) 


0.980 ± 0.014 


average (k) 


0.913 ± 0.012 



Table 7 

Ratios of the hadronic cross-section ratios at 137 GeV/c and the nuclear cross- 
section ratios around 600 GeV/c. 

where k is a parameter specific for the cross section ratio (compare with ta- 
ble 7). If we set K = 1 for simplicity, we see that the precision of (37) is about 
10%. The precision is improved by adequate adjusting of k. 

Unfortunately we cannot empirically derive k from experimental cross sec- 
tions for laboratory momenta around 600 GeV/c as necessary cross-section 
data is missing. Thus, as we want to deduce hadron-nucleon cross sections 
from nuclear cross-section ratios with best precision, we improve the rela- 
tion between hadron-nucleon and hadron-nucleus cross sections using the total 
cross-section model that was introduced in section 7. 

The idea of the model-based ratio method is the following: Rewriting (37) 
yields the following relation between the experimental hadron-nucleus and the 
model based hadron-nucleus cross-section ratios. 



^tot(pA) 



a. 



mod 



(A,atot(hN)) 



--mod 
Otot 



(A,atot(pN)) 



(3J 



experimental theory + (Jtot-data 



Taking the ratio of model based quantities reduces the effect of uncertainties 
in the cross-section model. Further, as data for o'tot(pp) is most precise and 
exists over a large energy range, it is convenient to use proton-nucleus cross 
sections in the denominator. As the energy dependence of the pp cross section 
is known at SELEX energies and the model is adjusted to describe NA cross 
sections for piab > 10 GeV/c, the energy dependence of (Ttot(hp), which we 
want to determine, is taken into account through the energy dependence of 
the pp cross section. 

To deduce the cross section (Ttot(hN) from the measured nuclear cross-section 
ratio, we fix (Ttot(pN) (= o'tot(pp)) first and calculate the denominator 
o"^°'^(A, (Ttot(pp)) by taking o'tot(pp) from parameterization (27) evaluated at 
the laboratory momentum of the nuclear cross-section ratio as given in table 6. 
Iterating with respect to the model input parameter (Ttot(hN) until the model 
based total cross-section ratio in (38) equals the experimental one, yields the 
desired hadron-nucleon cross section. At SELEX energy we interpret the result 
crtot(hN) identical to crtot(hp). 
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8.2.1 Results for atoti^ N) and atoti'^ N) using the ratio method 

Results of the ratio method are presented in table 9 together with the re- 
sults from the CH2 - C method. The errors of hadron-nucleon cross sections 
resulting from the ratio method include both the error in the measured nuclear 
cross-section ratio and model uncertainties. Model uncertainties are taken into 
account by adding the error of a model cross-section ratio in quadrature to the 
error of the corresponding experimental cross-section ratio given in table 6. 
The error in the model cross-section ratio is derived from the discrepancy be- 
tween model and measured cross sections observed for pA and vr^A total cross 
sections. Typical sizes of these discrepancies are shown in table 8. 





measured cross 


calculated 


cross-section 


nominal 


reaction 


section x fcgcaic 


cross section 


difference 


Plab 




[mbarn] 


[mbarn] 


[mbarn] 


[GeV/c] 


crtot(7r"Be) 


188.7 


188.8 


0.1 


640 


<7tot(vr"C) 


234.1 


231.4 


2.7 


590 


atot(pBe) 


270.2 


277.0 


6.8 


640 


o-tot(pC) 


336.8 


335.9 


0.9 


590 



Table 8 

Discrepancy between model and measured total cross sections. The measured pA 

cross sections are scaled by /cscaio- 

Further, as two different parameterizations for (d^o"/dtdM^)|t=o are available, 
we evaluate the ratio method for both, average the results and include their 
difference in the error of the mean. 

Finally, we want to mention that as little data exists for S^ scattering, we 
insert in the computation of (j™°'^(S^A) for B^-n, p's-n ^^^ (d^cr/dtdM^)|t=o, 
the parameterizations from pp-reactions. 



method 
description 


cTtot(S-N) 
[mbarn] 


cTtot(vr-N) 
[mbarn] 


Plab 

[GeV/c] 


difference method 


33.7 ± 3.1 


26.0 ± 2.1 


585 


ratio method. Be data 


37.4 ± 1.3 


27.1 ± 1.5 


640 


ratio method, C data 


37.0 ± 0.8 


26.4 ± 1.3 


595 


total average 


37.0 ± 0.7 


26.6 ± 0.9 


610 



Table 9 

The total cross sections (Ttot(5^ N) and (Ttot(7r^N) resulting from all methods and 

their average. 

Comparing the hadron-nucleon cross sections of the ratio and the difference 
method, we find good agreement of the results with respect to their errors. 
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As final result, we average the hadron-nucleon cross-section values from all 
methods. These total averages are presented in the last row of table 9 together 
with a corresponding averaged laboratory momentum. 



8.3 Comparison to models 



8.3.1 Comparisons for atotij^ p) 

Most of the models and parameterizations for hadron-nucleon cross sections 
exploit the interplay of 2 contributions. The Pomeron contribution, which 
dominates asymptotics at high energies; and the Regge contribution, which is 
important at low and medium energies. Many models (e.g. [35,36]) describe 
the energy dependence of total cross sections quite well. Because of this, we 
simply display in figure 14 experimental data from [8] and SELEX together 
with the parameterization for o"tot(7r"p, s): 



h0.079 



atot(vr"p,s)=35.9s"'^-*^ + 13. 7s^ 

for piab > 10 GeV/c, o"tot in mbarn, s in GeV^ , 



(39) 



which was presented in the particle data book 1996 [37]. We want to mention 
that in the data files of the particle data group errors in the high precision 
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fit of PDG 1996 (39) 
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Fig. 14. Existing data for (Ttot(7i" p) in comparison with our results and parameter- 
ization (39) of the particle data group 1996. 

data of A.S. CaroU et al. [38] are not completely taken over from the reference. 
Thus, we added the missing error contributions (systematic error due to target 
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density and extrapolation uncertainty) in quadrature. 

We point out that so far the total cross section o"tot ('?'""?) has been measured 
only up to piab = 370 GeV/c [38]. Thus, the SELEX total average for o"tot(7i"~N) 
at 610GeV/c is the first new measurement at higher laboratory momentum. 

In figure 14 parameterization (39) of the Particle Data Group, which uses a 
Pomeron intercept of 0.079, is overlaid to the data. The qualitative inspection 
of (39) suggests that it is strongly weighted by the huge amount of low energy 
data points and does not sufficiently well take into account the very accurate 
data of [38] at high energy. Our result seems to strengthen the trend observed 
in data of [38]. This trend implies a faster rise of the vr^p cross section with 
increasing energy than represented by (39). We just want to point out this 
observation, which may turn out to be in conflict with the belief that the 
energy increase of hadronic cross sections is universal. Further, we do not give 
any quantitative estimate of the Pomeron intercept for the 7r~p cross section. 
Its value is correlated to the assumed Regge contribution at low energy and 
its determination requires a careful analysis of the data. 



8.3.2 Comparisons for atoti^ p) 



Up to now, data on the total cross section crtot(S~p) are still scarce. In 
the past, there have been only two hyperon-beam experiments [4,39] giving 
information about the behavior of o"tot(S~p) in the momentum range 19 - 
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Fig. 15. Existing data for o"tot(S p) in comparison with our results and predic- 
tion (40). 

136.9 GeV/c. The SELEX total average for o'tot(S^N) provides first new data 
at higher energy. Figure 15 shows a compilation of data from previous experi- 
ments together with the SELEX result. Our measurement is 2.9mbarn larger 
than the data point at 136.9 GeV/c from [4]. It shows the rise of o"tot(S~p) 
with increasing beam energy. 
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Overlaid to the experiment data is the prediction for atoti^ P, Piab) from H. 
Lipkin (see [36]), which is given by: 



atot(S-p,Piab) = 19.5 (^)°''' + 13.2 (^)'°'' (40) 

for piab > 10 GeV/c, atot in mbarn, piab in GeV/c . 

The corresponding curve in figure 15 shows good agreement between our mea- 
surement and this prediction. 

It would be certainly desirable to find the Pomeron intercept for the S~p 
cross section. But the lack of low energy data does not allow any reasonable 
estimate of the intercept. 



9 Conclusions 



The SELEX collaboration has measured the total cross sections cr^^ti'^ Be), 
o"tot(7r"C), o-tot(vr~Cu), o-tot(vr~CH2), o-tot(S~Be), o-tot(S~C), o-tot(S"Cu), 
o-tot(pBe) and (TtotlpC) in a broad momentum range around 600 GeV/c us- 
ing a transmission method that was adapted to the specifics of the SELEX 
spectrometer. The accuracy of the results is within 0.6% - 1.5% for Be, C and 
CH2 and about 17.5% for Cu. 

The ratios of hadron-nucleus cross sections for Be and C show that 7r~- 
nucleus cross sections are a about factor of 0.7 lower than corresponding 
proton-nucleus cross sections. Furthermore, we find that the S~-nucleus cross 
sections are about a factor of 0.92 smaller than corresponding proton-nucleus 
cross sections. 

We observe that the results for (Ttot(pBe), (Ttot(pC), (Ttot(7r"Be), (JtotlTr^C) 
and (Ttot(7r^Cu) join smoothly corresponding cross-section data at lower en- 
ergy. The good agreement of the proton-nucleus and the 7r~-nucleus cross 
sections to Glauber model calculations which include an inelastic screening 
correction and one adjustable parameter in the density distribution, justify 
the deduction of atoti^p) and o-totl^r^p) from the nuclear cross sections. 

We deduced the hadron-nucleon cross sections atotlTr^N) and o'tot(S~N), 
which we regard as o"tot(7r~p) and o"tot(S~p), from our nuclear data using a 
CH2 - C difference and a model based ratio method. Results from the dif- 
ference method have an accuracy of 8.1-9.2%, while results from the ratio 
method have an accuracy of 2.2-5.5%. 

The total averages of all methods represent first new measurements for 
o"tot(7r~p) and (Ttot(S^p) at higher energy. Our result for o-tot(S^p) shows 
clearly a rise of this cross section with increasing beam energy, which agrees 
with the prediction of [36]. 

Our result for crtot(7r^p) joins nicely the high energy data of [38]. As men- 
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tioned in section 8.3.1, the data of [38] and our result may indicate a faster 
increase of the 7r~p cross section than predicted by the parameterization given 
by the Particle Data Group in 1996. 

An indication of a faster increase of the 7r~p cross section compared to the 
pp (and pp) one should be verified by a high statistic measurement using a vr^ 
beam and a hydrogen target to avoid systematic errors inherent to the method 
used in this experiment. In our opinion a measurement of the 7r~p cross section 
at 600 GeV/c or higher is the only experimentally accessible case to test if the 
high energy behaviour of the hadronic cross section can be different from the 
pp and (and pp) one. 
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